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Abstract 

A method to construct trihamiltonian extensions of a separable system is 
presented. The procedure is tested for systems, with a natural Hamiltonian, 
separable in classical sense in one of the four orthogonal separable coordinate 
systems of the Euclidean plane, and some explicit examples are constructed. 
Finally a conjecture on possible generalizations to other classes of systems 
is discussed: in particular, the method can be easily adapted to the eleven 
orthogonal separable coordinate sets of the Euclidean three-space. 

1 Introduction 

Separation of variables for Hamilton-Jacobi equation is a very effective method 
to find solutions of Hamiltonian systems. A set of coordinates in which a 
certain Hamilton-Jacobi equation can be separated is called a separable set 
of coordinates for that Hamiltonian system. The classical characterization of 
separability is restricted to the particular class of Hamiltonians defined on the 
cotangent bundle of a Riemannian manifold (Q,g) with the form: 

H = l -9 lJ PlP3 + V(q) = G(q,p) + V(q) . 

This kind of Hamiltonians is usually called natural. Moreover, in the classical 
theory, the separation is always performed through a contact transformation 
of coordinates, i.e. a fibred symplectic transformation on T*Q obtained from 
a transformation in the configuration space Q. In this contest a natural 
Hamiltonian is separable in a given set of coordinates only if its geodesic part G 
is separable J2| ■ It is therefore possible to determine the sets of coordinates in 
which a certain Hamiltonian may be separated through the study of all separable 
coordinate systems associated to the metric g. The particular case in which the 
separation is performed in an orthogonal set of coordinates is called orthogonal 
separation, other cases are said of general separation. In recent years a rich 
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and coordinates-indipendent characterization of both orthogonal and general 
separation was developed (see ^3 El El an d references therein). 

The case of orthogonal separation is quite important not only because it is 
simpler, but also because it has been proved that in a Riemannian manifold of 
constant curvature only orthogonal separation is possible (see ^2] f° r a review of 
such classical results) . The classical analysis of orthogonally separable systems is 
due to Stackel and Eisenhart £3 . The central key of their analysis is the concept 
of Stackel matrix, i.e. a n x n invertible matrix S (where n is the dimension 
of the Riemannian manifold) such that each element Sf belonging to the i-th. 
row of the matrix depends only on the coordinate qi. Stackel proved that a 
geodetic Hamiltonian is orthogonally separable if and only if the nonvanishing 
controvariant components of the metric tensor, g t% , form a row of the inverse of 
a Stackel matrix: 

g u = (S- 1 )) for a fixed j. 

The existence of a Stackel matrix is equivalent to the existence of a family of 
n—l particular Killing tensors Kn\ for the metric g. The vector space spanned 
by this family and the metric is the Killing-Stdckel algebra associated to the 
orthogonal coordinate systems [3] . Vice versa, from a Killing-Stackel algebra it 
is possible to reconstruct the Stackel matrix of the coordinate system. Finally 
also the separability of the complete natural Hamiltonians can be related to the 
Stackel matrix: the Killing-Stackel algebra allows to construct a family oi n—l 
potentials such that the Hamiltonians 

H (i) = K (i)PiPi + v (i) 

are in involution both mutually and with H. An alternative way to restate these 
results is the following: a family of n Hamiltonians Hi is orthogonally separable 
if and only if there exists a Stackel matrix S and a set of n functions $>i(qi,Pi), 
each depending only on the i-th coordinate, such that 



S 









H 



$i(<Zi,Pi) 



These classical results has been put in a more handy form in pQ, where it 
has been proved that a natural Hamiltonian is orthogonally separable if and 
only if exists a Killing tensor K with n normal and simple eigenvalues such 
that d(KdV) = 0. A stronger result can be proved if Q is also equipped with 
a conformal Killing tensor with vanishing Nijenhuis torsion, i.e. a tensor L of 
type (1,1) satisfying 

[LX, LY] - L[LX, Y] - L[X, LY] + L 2 [X, Y] = 

for all vector fields X and Y on Q, and 

{L^pip^g^piPj} = cg lJ piPj 
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In this case it is indeed possible to generate the whole Killing-Stackel algebra 
by a recursion formula. These systems are called Benenti systems or L-systems. 
The recurrence property of Benenti systems was explained in the paper [5], 
where an extension of the systems with a Killing-Stackel algebra generated 
by the conformal Killing tensor L is constructed, and it is proved that this 
extension fits in a bihamiltonian hierarchy defined on an extended manifold 
that contains the original phase space T*Q as a symplectic leaf of a degenerate 
Poisson structure. Therefore the recursion formula for the Killing tensor in the 
Killing-Stackel algebra is obtained from the bihamiltonian recursion scheme. 

So far the only contact transformations has been considered in order to 
perform a separation of variables. In 8 Sklyanin proposed to generalize this 
framework allowing general symplectic transformations and substituting the 
separation relations expressed by Q with the n general equations: 



Also this more general kind of separation has been interpreted in the bihamil- 
tonian framework [5]. A different point of view could be found in the review 



The aim of the present work is to introduce, in the simple case of the Eu- 
clidean plane, a recipe for constructing trihamiltonian (and consequently bi- 
hamiltonian) extensions of a orthogonally separable natural Hamiltonian. This 
procedure works also for systems separable in symmetric sets of coordinates, 
whose Killing-Stackel algebra is in general not constructible from a conformal 
Killing tensor L. Although the method seems promising also for the separable 
coordinate systems of the Euclidean three-space, the possibility of its extension 
to any orthogonally separable system remains unfortunately to be proved in 
general. 

The construction presented in this paper in not interesting just as a new 
class example of trihamiltonian systems, in which the third Poisson structure 
links together different Lenard chains. Indeed it gives also an Hamiltonian 
interpretation of the concept of separation curves presented in [2]: this 
curve is obtained through a reduction procedure from the common Casimir 
function associated to the trihamiltonian structure. Moreover starting from 
a trihamiltonian extension of a separable system in the plane it is possible 
to construct two different bihamiltonian extensions, whose recursion relations 
are in a case of "unsplit" and in the other of "split" type. This suggests the 
possibility to deal with both cases in a unified way. 

The paper is organized in the following way: in the next section the four 
separable set of coordinates on the Euclidean plane are briefly presented, as well 
as the pair of quadratic Hamiltonians associated with each coordinate systems 
and the corresponding Stackel matrix. Because the Euclidean plane is a constant 
curvature manifold, separation of variables only occurs in orthogonal coordinates 
and therefore the presented list is exhaustive. 
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Section three first presents the procedure for constructing the extended 
systems in the two asymmetric sets of coordinates, moreover some basic concepts 
of the bihamiltonian and trihamiltonian framework are reviewed. In the 
following the procedure is separately adapted to each of the two symmetric 
coordinate systems; this allows to show how this method could be generalized 
to other cases. 

Finally, some explicit examples are shown. In particular the one-Casimir 
extensions of the Henon-Heiles and Kepler systems presented in [2] are recovered 
from a trihamiltonian point of view and it is shown, in the Kepler case, how 
multiple sets of separation coordinates lead to different trihamiltonian extension 
of the same Hamiltonian. 



2 Systems related to separable orthogonal web 
in the plane 

It is well known that, in the Euclidean plane, there are only four systems 
of orthogonal coordinates allowing to separate the Hamilton- Jacobi equation 
associated to a natural Hamiltonian As a matter of fact, instead of 

considering the coordinate system, it is often appropriate to consider the 
associated web, i.e. the family of curves on which the coordinates are constant 
because the web associated to a coordinate systems is invariant respect to 
transformation qi ^ Qi such that the new coordinate Qi turns out to depend 
only on qi for all fixed i. This kind of transformation of coordinates is called 
separated. If x and y denote the Cartesian coordinates in the Euclidean plane, 
then the other three coordinate systems are: 

• parabolic coordinates: 



x = u = x + yjx 1 + y 2 
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y = \/—uv v — x — \/ x 2 + y 

whose associated web is made up of confocal parabolae with focus in the 
origin and symmetric respect to the x-axis; 

• elliptic- hyperbolic coordinates: 



2A- 



z=§ s = ^{x + k/2) 2 + y 2 + y/(x - k/2) 2 + y 2 

y= VHfEgEEn d =^/(x + k/2) 2 + y 2 - ^{x - k/2) 2 + y 2 

whose web is made up of a family of confocal hyperbolae and a 
family of confocal ellipses, whose common foci lie on the x-axis with 
coordinate ±fc/2; 

polar coordinates: 



x = r cos 6 r = \J x 2 + y 2 
y = r sin 9 9 = arctan(y/ir) 
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centred at the origin. 



The study of the Killing tensors associated to each coordinate system allows 
to construct four canonical types of completely integrable systems (see [3] and 
references therein). In fact, a pair of Hamiltonian functions H and K, in 
involution with respect to the canonical Poisson bracket on the cotangent bundle 
of M 2 , and both separable in the given coordinates, is related to each coordinate 
system (more precisely to each web underlying the coordinate system). The 
separability of these Hamiltonians is ensured by the following relations involving 
the Stackel matrix of the coordinate system. 

Cartesian coordinates 

Hamiltonian functions: 

H = -pl + Mx) 

k = -P 2 y + Mv)- 

Stackel relation: 





Parabolic coordinates 

Hamiltonian functions: 

H = 

K = 
Stackel relation: 



Elliptic-hyperbolic coordinates 

Hamiltonian functions: 



[\pI + M u )) u - {\pI + <fo( v )) v 





H 

K = 
Stackel relation 



M + ^fl) ( s2 ~ fc2 ) ~ M + ^ {d) ) {d2 - fc2) 

s 2 -d 2 

{hp* - M + tew - frtfS ^ - fc2 ^ 2 - fc2 ) 

s 2 -d 2 




±p 2 d + Md) 
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Polar coordinates 

Hamiltonian functions: 

K = -\p 2 e -fc{8). 
Stackel relation: 

1 ^ \( H ) = ( ^ Pl + Mr) ) 

o 1 J \ K J \ -hti-M*) J' 

Before going ahead a remark on the choice of the Hamiltonians is needed: 
sometimes, in Cartesian coordinates, the two Hamiltonian functions 

R = -(pi + P 2 y ) + + My) 
k = \p 2 y + Mv) 

are considered, with the Stackel matrix 




The Hamiltonians H and K used in the present work are trivial recombi- 
nations of Hamiltonians H and K that allow to put the Stackel matrix in a 
more useful form. This form is such that in the i-th row there is just an ordered 
sequence (possibly decreasing) of powers of a suitable function of the coordinate 
qi. One of the key point in order to extend the procedure presented in this paper 
to more general cases is to characterize the class of Stackel matrix that can be 
put in this form. 

It is worth to observe that, with an appropriate change of variables, the 
previous four types of separable systems are mapped into the four types of 
integrable systems in the plane with two quadratic first integrals (see, for 
example 0). The method to construct a trihamiltonian extension of the 
previous systems is slightly different for the asymmetric coordinate systems 
(parabolic and elliptic-hyperbolic) and for the symmetric ones (polar and 
Cartesian). This difference is due to the different form of the Stackel matrix 
in the two cases, as it can be seen from the above formulae. In fact, for 
asymmetric coordinate systems all the rows of the Stackel matrix have the same 
form and differ only in the current coordinate, whereas for symmetric coordinate 
systems different rows of the Stackel matrix have different forms. Therefore, the 
method that will be used in the asymmetric case will need to be adapted to each 
symmetric coordinate systems in a specific way. 
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3 Trihamiltonian extension of separable systems 
3.1 Asymmetric coordinate systems 

The idea of constructing a trihamiltonian extension of separable systems 
originates from the construction of bihamiltonian extension of Benenti systems 
presented in [BJ, and it is motivated by some results about trihamiltonian 
systems presented in |3] that will be briefly summarized. If P, Q and R are three 
Poisson structures mutually compatible on a manifold of dimension 2n + k, and 
/ is simultaneously a Casimir function for the two Poisson pencils Q — XP and 
R — fiP, polynomial in A and /i, then all the coefficients h% of this polynomial arc 
functions mutually in involution with respect to all the three Poisson structures, 
and satisfy a recurrence scheme depending on the form of the polynomial /. 
Under some further hypotheses, it is possible to find out a set of coordinates 
{Xi, [ii,c a } with the following characteristics: 

1. c a are k Casimir functions for P 

2. Xi, Hi are 2n Darboux-Nijenhuis coordinates for both the Nijenhuis tensors 
obtained, thanks to a deformation procedure, from the two Poisson pencils 
Q- XP e R- fiP 

3. Xi,/ii satisfy the n relations: 

f(\um; {h k }, c a ) = n(X u iJ,i) (2) 
in which the functions n(X, fj,) are constant functions on the manifold. 

The relations J5J imply the separability (in Sklyanin sense) of Hamiltonians hk ■ 
It is then natural to ask if it is possible, starting from a classical separable 
system, to obtain an adapted trihamiltonian structure. 

As a matter of fact, the separation relations encoded by the Stackel matrix 
in asymmetric coordinate systems are very similar to relations , in the special 
case in which the dimension of the symplectic leaves of P is n = 2. Indeed with 
the symplectic (and separated) transformations: 

Ai = u, A 2 =v,fii= p u , M2 = Pv 
for parabolic coordinates and 

Xi = s 2 — k 2 , A 2 = d 2 - k 2 , Hi = ^, fi 2 = 7n 

Zs Id 

for elliptic-hyperbolic coordinates, the relations extracted from Stackel matrix 
become: 

/Ax 1 \( H \ ( ri(A 1)Ml ) \ 
\X 2 I )\K ) \ r 2 (Ai, M i) J 

that is 

K + X l H = r l {X t ,^) (3) 
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where the functions ri(A,/i) are, respectively in the two coordinate systems: 



t(A,/x) = -A^ 2 + A0 2 (A) 



r,(A,/i) = 2 A i 2 (A 2 + Afc 2 ) + A4(v / A + fc 2 ). 

The difference between relations (0 and is the absence in the second ones 
of any Casimir function c a ; this implies that the polynomial K + XH can not 
be a Casimir function for any Poisson pencil, including the canonical Poisson 
structure on the cotangent bundle of M 2 . In order to obtain a set of relations 
analogous to J5J it is necessary to extend the phase space of the system with 
a suitable number of coordinates c Q , and to define on the extended space a 
Poisson structure P such that the new coordinates are its Casimir functions, 
and the reduction of P on its symplectic leaf {c a = 0} gives the original Poisson 
bracket. This allows to think relations 10 as a reduction of relations © on 
the symplectic leaf {c a = 0}. Moreover, in the extended space, the functions 
that reduce to H and K could depend on the new coordinates c Q , hence it is 
also necessary to introduce two "deformed" Hamiltonians H and K. Taking in 
account these observations the simplest function analogous to the one appearing 
in J2J and compatible with the relations @ is: 



/ = K + XH + ciA 2 + /ic 2 + \[ic 3 



(4) 



where H and K are the appropriate extensions of the Hamiltonians H and K 
and three extra coordinates are introduced. The extended Poisson structure P 
has hence the form 

/0 10000\ 
1 
-1 
-10 


\ 00000/ 



p 



And the trihamiltonian recursion scheme associated to the function is: 



C2 c 3 



K 




In the previous scheme the notation / — > Xf means that the vector field Xf is 
obtained from the differential of the function / through the Poisson structure P. 
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Further, if the relation Pdf = Qdg between the two functions / and g holds, it 
is sometimes said (understanding some privileged role for the Poisson structure 
P) that / is the bihamiltonian antecedent of / through the structure Q. 

The previous construction needs some remarks: first, it is possible to reduce 
the number of new coordinates that is needed by omitting the "vertical" part 
of the recurrence scheme and looking for a function in the form: 

f = K + XH + A 2 ci . 

This approach is followed in |U| and [2], but allows to obtain just a bihamiltonian 
recursion. On the other hand it is possible, both in the bihamiltonian (as in [2]) 
and in the trihamiltonian case, to add to the function Q a "redundant" Casimir 
function, that is a Casimir function common to all the Poisson structures, using 
the function: 

/ = / + A"c 4 with n > 2 . 

In this way a different pair of extended Hamiltonians H and K is constructed 
and it can be reduced to that previously obtained putting C4 = 0. Even the new 
Poisson structures can be related to the old ones through the reduction on the 
zero level set of the common Casimir function 04. 

Besides these digressions, the next step of the procedure is to require that 
the coordinates A^ , /ii are separation coordinates in the sense of (0) . This implies 
that H and K must solve the linear system 

f(\i,Hi;K,H,c a ) = n(\i,fj,i) 

Therefore the solution, using the function is: 

jj = r 2 (\ 2 ,n2) - nXWn) _ 

A2 — Ai 

(Ai + A 2 )ci + — — c 2 H c 3 ; 

Ai — A2 Al — A2 

A 2 ri(Ai,/xi) - Air 2 (A 2 ,/i2) 
A2 — Ai 

, , . (A 2 mi-^2Ai) A 2 Ai (fi 2 - 

A2A1C1 H C 2 H r C3 . 

Ai — A2 A2 — Ai 

The last step, in order to obtain the desired recursion scheme, is to construct 
other two Poisson structures, after the P defined above. Following 0] one start 
defining the two "deformed" Poisson tensors represented by the 7x7 matrices 



Qd - 



K 



Rd — 






A 





-A 




















M 




—M 
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where 



A = 
M = 

and then the two vector fields by constructing 



Ai 

A 2 

Mi o 




x q - J2 a F <*-£z 

where the functions F a and G a are the bihamiltonian antecedents of the Casimir 
functions of P through, respectively the tensors Q e R. This means that 
relatively to the recursion scheme (JjjJ, they satisfy the relations PdF a = Qdc a 
e PAG a = Rdc a . Hence in the given case it holds F-y = H, F2 = F3 = and 
Gi = 0, F 2 = K, F 3 = H . The two vector fields then are: 

X R = K— + H-. 
Lastly, it is possible to construct the two tensors 

Q = Q d -L XQ P 
R = Rd — Ljfjj P 

The surprising result of this complicated construction is the following fact, that 
could be verified by direct computation: 

Fact 1 The tensors Q and R previously constructed are compatible Poisson 
tensors independently of the choice of the functions rt(X, fi). Moreover, the 
function f given by ^j) is a common Casimir function for the two Poisson 
pencils Q — XP and R — XP, realizing the recursion scheme $Bjj. 

Being ci and C3 two common Casimir functions of P and Q, both these 
Poisson structure are reducible on the level surface {02 = C3 = 0}. On this 
surface the function @ reduce itself to 

K + XH + A 2 ci 

and the trihamiltonian recursion scheme 10) becomes a simply bihamiltonian 
one. Because the separation relation (j2J for the reduced systems becomes 

K + XiH + X\ci = n(Ai,Mi) 
k + X 2 H + X\c 1 = r 2 (A 2 ,M2) 
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when <pi — cf>2, and then r\ = r 2 = r(X,n), the function 

K + XH + A 2 ci - r(A, fi) 

define an unsplitted separation curve in the sense of 

On the contrary, reducing the function / and the Poisson structures P and 
R on the level surface {c\ = 0} of their common Casimir function c\, a splitted 
bihamiltonian recursion chain is obtained. 



3.2 Symmetric coordinate systems 

The procedure described above cannot be directly applied to Hamiltonians 
separable in a symmetric coordinate system, i.e. in polar or Cartesian 
coordinates. In these cases, actually, the relationships obtained from the Stackel 
matrix cannot be summarized with an unique polynomial, able to suggest a 
candidate / to the role of common Casimir function to the two Poisson pencils 
Q — XP and R — XP. This difficulty could be overcame observing that a 
trihamiltonian structure doesn't admit necessarily only one common Casimir 
function for the two Poisson pencils. Therefore, it will be sufficient to build a 
different common Casimir function for each of the different polynomial relations 
generated by the Stackel matrix. 

In the case of polar coordinates, through the symplectic transformations: 

Ai = r 2 , X 2 = tan#, /ii = fi 2 — cos 2 9pg 

2r 

the Stackel relationships could be rewritten with the two equations: 

K + XiH = ri(Ai,/Xi) (f .s 
K - r 2 (A 2 , M2 ) W 

where the functions fj(A, fx) are: 

n(A,/i) = 2^ 2 A 2 + A0i(VX) 

r 2 (A,/i) = -i(l + A 2 )V -0 2 (arctan(A)) 

Hence two functions have to be constructed, with the form: 

/i = K + XH + A 2 ci + fic 2 + X^c 3 . , 

h = K + X Ci + fic 2 1 ' 

and realizing the trihamiltonian recursion scheme: 
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The expressions for H and K can be obtained by solving the system: 

\ K+ X 1 H + \\ci + ^ic 2 + Ai^ic 3 = ri(Ai,/ii) 
| K + A 2 c 4 + Ai 2 c 2 = r 2 (A 2 ,/i 2 ) 

whose solution is: 

~ ri(Ai,/xi) -r 2 (A 2 ,M2) . . M2 - Mi , A 2 

.H = A1C1 H c 2 - ^ic 3 + — c 4 

Ai Ai Ai 

if = r 2 (A 2 , M2) - M2C2 - A 2 c 4 

The P, Q and i? structures are obtained, analogously to the asymmetric case, 
by adding to the tensors Qd and Rd the Lie derivative of P respect to the 
vector fields Xq and Xr, constructed with the antecedents of the four Casimir 
functions of P. From the recursion scheme |JSJ these two vector fields turn out 
to be: 

Xq = H— + K— 

OCi OCi 

oc 2 dc 3 

By direct calculation the following fact can be proved: 

Fact 2 The tensors Q and R previously constructed are compatible Poisson 
tensors, independently of the choice of the functions rj(A, /i). Moreover, the 
functions f\ and / 2 given by ^ are two common Casimir functions for the two 
Poisson pencils Q — XP and R — XP, thus realizing the recursion scheme pj) . 

Eventually, in the case of Cartesian coordinates, the starting point is the 
construction of two functions of the form: 



/1 = H + CiA + pLc 2 
h = K + Ac 3 + LIC4 

that realize the trihamiltonian recursion scheme: 
c 2 c 4 



(9) 



X H X K 

H c i K c 3 



(10) 



By solving the linear system: 



H + C1X1 + [i\c 2 = ri(Ai,/ii) 
K + A 2 c 3 + ^ 2 c 4 = r 2 (A 2 ,^ 2 ) 
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it is possible to find out the Hamiltonians: 



H = ri(Ai,//i) - ciAi - /iic 2 
K = r 2 (A 2 ,/i 2 ) - A 2 c 3 - /i 2 c 4 

and from the recursion scheme (|10fl the following vector fields are obtained: 

~ d ~ d 
X Q = H—+K— 
oci oc 3 

d d 
X R = H—+K-. 

ac 2 0C4 

Even in this case it holds: 

Fact 3 The tensors Q and R previously constructed are compatible Poisson 
tensors, independently of the choice of the functions n{X,fJb). Moreover, the 
functions f\ and / 2 given by are two common Casimir functions for the two 
Poisson pencils Q — XP and R — XP, thus realizing the recursion scheme \TQjj . 



4 Examples 

4.1 Trihamiltonian extension of Henon— Heiles system 

An application of the recipe presented above is given by the Henon-Heiles 
system, already considered in 0, whose Hamiltonian expressed in natural 
coordinates {qi, q 2 ,pi,p 2 } is: 

H =\{p\+pl) + \qiql + ql 

This system is separable in parabolic coordinates and the previous Hamiltonian 
can be obtained from the general expression through the coordinate transfor- 
mation: 

qi = Ai + A 2 , q 2 = 2V-AiA 2 

pi - a!-a 2 ' - v-xiMf^x; 

and choosing the two arbitrary functions <f>\ e <j) 2 to be: 

cf>i(z) = 4> 2 {z) = z 3 
Applying the explained procedure one get the two deformed Hamiltonians 

H = -(pj + pt) + -qiql + qf - qict - 2—c 2 - ptc 3 

Z Z q 2 

K = o(<?2PiP2 - qiP 2 ) + T<7i<7 2 + tx<7 2 - -A%c\ + (2 Pi)c 2 - -<M5 2 c 3 

2 4 lb 4 q 2 Z 
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that are put in the recursion scheme 10 by the three Poisson structures: 
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where X = = PdH and Y = X^ = PdK. It is worth to observe that both 
P and Q can be reduced by restriction on the level surface C2 = C3 = and 
that the reduced Hamiltonians and Poisson structures are the ones considered 
in Instead the structure R isn't reducible, so the reduced system appear to 
be just a one-Casimir bihamiltonian extension of the Henon-Heiles system. An 
advantage of the knowledge of a trihamiltonian extension of the system is the 
possibility to perform a reduction on the level set c\ = 0, producing a different, 
two-Casimir bihamiltonian extension of the system, whose Poisson structures 
are the reduction of P and R. 



4.2 Trihamiltonian extension of Kepler system in the 
plane 

The Kepler system in the plane is separable in three different set of coordinates: 
the parabolic, the elliptic-hyperbolic and the polar coordinate systems. To each 
of these sets of coordinates can be associated a different second constant of 
motion (related to the Killing tensor of the coordinate system) and a different 
trihamiltonian extension. 

In the case of parabolic coordinates the transformation from the natural 
coordinates to the coordinates {Xi, Hi] is given by: 

<7i = 2 V-Ai A 2 , <? 2 = Ai + A 2 

„ _ v / -AiA 2 (jJi-jJ2) _ Aijji-A 2 M 2 

r 1 ~ A1-A2 ' V2 — A1-A2 
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and setting 



M*) = 



a 

2? 



a 



the deformed Hamiltonians H and K become 

a 



H 



K 
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They are put in the recursion scheme (JSJ by the three Poisson structures 
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reducing the Poisson structures P and Q and the Hamiltonians on the level 
set C2 = C3 = 0, it is obtained the one-Casimir extension considered in But 
it is important to observe that the separation coordinates used don't seem to 
originate a unique "separation curve" because the two functions <pi and <f>2 are 
different. 

In the case of the elliptic-hyperbolic coordinates the coordinates transfor- 
mation from the natural coordinates is: 



92 - k — 



Pi 



fe(Ai-A 2 ) ' 



_ V(Ai+fc2)(A 2 +fc2) 

9i — fc 



P2 



o Mi(Ai+fc )-M2(A 2 + fc ) / \ — \- 

2 fe(Ax-A 2 ) V-A1A2 
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and the two arbitrary functions are 

<M*) = <fo{ z ) 



k 2 



The two deformed Hamiltonians H and K turn out to be: 



H 



K = 



1 , 2 2\ 11 



1 



(q( + qi - 2 qi k) Cl + ^ 
1 



Pi 



(qi - k)q 2 pip 2 
1 



+- 



2fc 2 V (qi - k) 

1 

" 2 
fc) 2 



, 

— C 2 

92 



"P2 - 



2A:)p; 
(ft 



2 2 
9 2 Pl 



Pi 



C3 



2(gi - fc) 
afcgi 2 2 



2fc 2 V 92 
and the three Poisson structures are: 



(91 - k) 



-Pi 



1 

C2+2 



g 2 fc ci 



gi - fc 



- 92P2 C 3 



/ 








1 











\ 
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-i 
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Q 



o o 



qi (2k — qi) (fc — <ji)<?2 

(k-qx)q 2 -q\ 

—Xi —X'2 






(k 



qi (qi - 2k) 
(qi - k) q 2 


— qi)p2 + 92P1 

-x 3 






(qi -k)q 2 Xi \ 

q\ X 2 

(qi — k)p 2 — q 2 p\ X A 

X 4 

-X4, 



J 



R = 






(g| + fc 2 )pi 
2k'- ! (k-q 1 ) 
92P2 



-91 )P2 + 92P1 
2k 2 




-Yi 



(fc — gi)pi , 



91 (91 — 2k )P2 
292 k' 2 



-Y 2 3 
-X 2 



(g|+fc 2 ) P1 
2k' 2 ( qi -k) 
_ 92 P2 



(g 1 -fc)p 2 -92Pl 
2p 



(91 -fc) 2 P2+Pl92 
2fc 2 <j 2 (91 -fc) 





-Y 3 
-x 3 



(91- 



fc)P2-92Pl 



(91 -fc)Pl 
Si- 2 



91 (91 ~ 2fc )P2 



(91- 



2,2 ft 2 

fc) 2 P2 + Pl92 
T, — 



292*= (91-fc) 
P2P1 





-Y 4 



Yi X! 



Y 2 X 2 



x 3 



Y 4 X 4 



where X 



PdH and Y 



PdK. 
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Finally in the case of polar coordinates the transformation between the 
natural coordinates and the coordinates {Ai,/ii} is given by: 



91 
Pi 



Ai 

2/JiAi-^2A 2 -A 2 jJ2 



VAi(l+Al) 

and with the choice 

& (*) 

the deformed Hamiltonians become 
1 



92 = A 2 

P2 



Ai 



1+Af 
/J2A2+2A1 A 2 /ji+/J2 
VAi(l+Al) 



A" 



= 



(<7i + Ql) c i 



3lPl + 92P2 - 2qf(qip 2 - q 2 pi) 



2(qf + qlY 



1 



(q 2 2 Pt + qfpi) 



-c 2 



2 91P2 

9i 



5iPi + 172P2 



2(9? 

- 92P1 



-ci 



q\Y 

?i 



c 3 + 



92 C4 



9l(?i + <?i) 



-c 4 . 



They are put in the recursion scheme © by the three Poisson structures 
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li+i'i 
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J^r X1OOY1 \ 
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-X A 
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X2OOY2 
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0000 
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where X = = PdH, Y = X R = PdK and 



A 



C 



0i 



2 qiPi + Q2P2 + iq\(q\V2 - qiPi) 



,2 



i{ql + qlY 

qiPi + q2V2 - 2q\{qiP2 - 32P1) 



/' = <l\<h .; 

- ' 2 , „2\ 2 



2(<? 2 + <z 2 2 ) 

gf(giPi + 02^2) + gg| (gip 2 - ggPi) 

+ 92P2 - 2q\(qip 2 - g2Pi) 



= (qiP2-q2Pi) — -2 

2 (9? + <z 2 2 ) 



5 Final remarks 



In this article a procedure to construct trihamiltonian extensions of classical 
separable systems has been presented. This procedure has been tested in the 
particularly simple case of the Euclidean plane, but it is in principle applicable 
to more general cases. It is articulated in the following steps: 

1. Write down the separation relations of a given system with n degree 
of freedom in the form involving the Stackel matrix of the separation 
coordinates: 



/ an(qi) 
021(92) 

V am(q n ) 



ain(qi) \ 

a2n(<72) 



( #i({pi, 9j }) \ 
K 2 ({p i ,q l }) 



(q n ) ) \ K n ({pi, qi }) ) 



( $iOi,gi) \ 
$202,92) 

V $n(Pn,qn) J 



2. Transform, through a transformation to a suitable system of coordinates 
{Xi,/ii}, and in some case a linear combination of the Hamiltonians, the 
Stackel relation in such a way that the i-th row contains just ordered 
powers of A^. This is the most troublesome step, in fact it's not clear which 
kind of Stackel matrix can be put in this form, although it is possible for 
all the eleven orthogonal separable webs in R 3 . 

3. Now the separation relation encoded by the i-th row of Stackel matrix is a 
polynomial in the coordinate Xi , its coefficients are a set of Hamiltonians 
Hj obtained from the recombination of the Kj. These polynomial can 
be grouped on the basis of their form: to each different form of the 
polynomials it is associated a trihamiltonian recursion scheme between 
the Hj and a corresponding polynomial function. In general there will 
be m < n of such functions, labelled /1, . . . , f m . The coefficients of the 
polynomial functions // are the deformed Hamiltonians Hj , together with 
a suitable number of Casimir functions. 
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4 



Solving the n equations (linear in the n deformed Hamiltonians) 



/i(Ai;{^}) 



< 



k fm{\i\{Hj}) 



the explicit forms of the deformed Hamiltonians is obtained. 

5. Extending in a trivial way the canonical Poisson tensor, the Poisson 
tensor P on the extended space is obtained. Instead, from the deformed 
Hamiltonian and the recursion scheme, it is possible to construct the two 
vector fields Xq and Xr, and then the two tensors Q = Qd — ^x Q ,p and 
R = Rd — Lx R ,p ■ 

6. The conjecture, on which this procedure is based, is that the three 
tensors P, Q and R constructed in this way are all Poisson tensors 
mutually compatible. Moreover the m functions /z(A, )i) are common 
Casimir functions for the two Poisson pencils Q — XP and R — [iP . 
As a consequence the deformed Hamiltonians Hj satisfy the required 
trihamiltonian recursion scheme. 
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